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ABSTRACT 

We report detection of strong circularly polarized emission from the transient 
bursting source GCRT J1745-3009 based on new analysis of 325 MHz GMRT 
observations conducted on 28 September 2003. We place 8 R as the upper 
limit on the size of the emission region. The implied high brightness temperature 
required for an object beyond 1 pc and the high fraction of circular polarization 
firmly establish the emission as coherent. Electron cyclotron or plasma emission 
from a highly subsolar magnetically dominated dwarf located <4 kpc away could 
have given rise to the GCRT radio emission. 

Subject headings: Galaxy: center - radio continuum: general - stars: individual 
(GCRT J1745-3009) - stars: variables: other 
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Introduction: 



The present generation of largely parabolic dish-based radio interferometers and data 
analysis procedures are poorly optimized for discovering transient emission. Consequently, 
the characteristics of transient radio source populations remain highly undersampled. 
Nevertheless, dedicated programs are revealing novel radio transients with unique emission 
characteristics, including GCRT J1745— 3009. Discovered in 325 MHz archival VLA data 
taken in September 2002, this source exhibited ~10 minute long, 1 Jy peaked bursts with a 



77 minute period over the ~7 hour baseline available in the discovery data (IHyman et al. 



20051 ). The measured time scales of the bursts implied a brightness temper ature exceedin g 



the 10 K Compton limit for a distance beyond ~100 pc. As discussed in 



Hyman et al. 



( 120051 ). its characteristics did not match any known mechanisms of emission in transient 
compact sources. As a result, GCRT J1745— 3009 appeared to represent a member of a new 
class of coherently emitting objects. 

Search to detect additional bursts from the source in archival data sets were 
made, resulting in its re-detection for ~2 minutes at 325 MHz Giant metrewave radio 



telescope (GMRT) o bservation in September 2003 ( IHyman et al. 



( IHyman et al. 



20061 ) and in March 2004 



20071 ). At the last known epoch of emission detected in 20 04, the source 



exhib ited an unusually steep spectrum with alpha 
20070 . 



•13±3 (S (i/)oc v a ) (IHvman et al. 



Several t h eories have been proposed to explain the emission from GCRT J1745-3009 (se e 



Hyman et al 



20051 ). double 



dwarf pulsar (IZhang fc Gil 



(200 7J) and reference s therein). T hese include mil" 



(ITurolla et a 



2005), precessing 



20051 ). Also, 



Zhu k Xu 



Hallinan et al. 



ing ( IKulkarni fc Phinney 



2006) and a transient white 



( 120071 ) suggested a nearby ultracool 



dwarf as its progenito r. Unfortunately, n on-detection of GCRT J1745— 3009 at frequencies 



other than 325 MHz ( iKaplan et al. 



20081 ) hindered attempts to discover its progenitor. Its 
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brightness temperature in possible excess of the Compton limit suggests coherent processes 
such as electron cyclotron maser emission or plasma emission. Notably, both of these 
mechanisms produce circularly polarized emission, motivating us to re-examine earlier 
detections of GCRT J1745— 3009 in search of its signature. We have reanalyzed the 325 
MHz GMRT data from September 2003, the detection with the highest signal to noise ratio 
per integration time available. The subsequent sections are arranged as follows - in Sect. 2 
and 3 we describe our reanalysis procedure and its results, respectively. The interpretation 
of these results are presented in Sect. 4 and conclusions in Sect. 5. 

2. Observations and data analysis: 

During the observations, the pointing center of the GMRT antennas was 0.5° West of 
GCRT J1745— 3009 (hereafter GCRT). The unpolarized calibrator 3C48 was observed for 
absolute gain calibration - see Hyman et al. 2006 for further details of the observation and 
total intensity calibration. At the time of these observations, the full polarization mode of 
the GMRT was unavailable and only parallel-hand correlations (Stokes RR and LL) were 
recorded. If the polarization leakage coefficients of antenna 'i' is denoted by D i? then the 
observed right circularly polarized signal on the i-j baseline is 

R[.R'* = (l-D i -D* j ).R i .IT j +R i .(D j L j )* + (D i A i ).ir j + 

(1) 

In general, the quantities are complex; we denote the 'true' values without primes and 
the measured quantities with primes; and a similar equation holds for the left circularly 
polarized signal L^L'*. During standard calibration without full polarization information, 
the leakage terms are considered small and discarded. However, the above equation shows 
that even for an unpolarized source the last two terms [(DiRi) .(DjRj)* + (DiLi).(DjLj)*] 
lead to amplitude closure errors. If the source is linearly polarized (Ri.L* ^ 0), additional 
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errors are added. 

To determine magnitudes of these leakage coefficients, the connections for the 2 
polarizations of a single antenna is reversed. This allows acquiring cross products of 
opposite handed polarizations of one antenna with the rest of the antennas. Observations 
of strong unpolarized calibrators have shown that D { < 0.1 at 325 MHz. Since the leakage 
terms for different antennas are expected to be uncorrected, the estimated error in Stokes 
V in the image plane for even a 100% linearly polarized source will be < 5% of Stokes I. 
Moreover, the polarization leakage in antennas cannot add coherently in this case, indicating 
absence of any significant spurious small diameter source in Stokes V in the image plane. 
However, the above analysis is valid for a source along the axis of the antennas and the 
GCRT was off-axis by 0.5°. Antenna off-axis positions may suffer from correlated patterns 
among the antennas (e.g., different FWHM of the right and left circularly polarized primary 
beams) leading to non Zero Stokes V. Therefore, we obtained supplementary polarization 
calibration observations of the unpolarized calibrator 3C48 at 325 MHz with 20 different 
angular offsets within the FWHM of the antenna primary beam. Fig. 1 shows that the 
pointing positions sampled a range of primary beam offsets extending beyond the GCRT 
offset as observed on 2003 Sep 28. These observations were sufficiently short such that 
rotation in feed parallactic angle within any particular pointing position were negligible. 
For all the off-axis positions after calibration, the maximum value of Stokes V was found to 
be 0.02 times Stokes I. Hence, the error in Stokes V within the FWHM of the antennas is 
<2% of Stokes I for an unpolarized source and <5% for a high linearly polarized source at 
325 MHz band of GMRT. 

A variant of multi-resolution Clean was used to better image the non-variable and 
extended field sources. After amplitude and phase self-calibration of non-variable sources, 
all sources except one were subtracted from the uv-data. and the inverse of the amplitude 



AZ 




Fig. 1. — 3C48 observed through 20 different parts of the primary beam. The outer circle 
indicates the FWHM of the primary beam (radius 45') and the inner circle is half of that 
distance. Azimuth is along X-axis, and elevation along Y-axis. The location of the GCRT is 
indicated by an open circle. 

self-calibration table applied to avoid changes in the flux density of the GCRT during 
imaging. Subsequently, the GCRT and the strong compact field source (G358. 917+0. 073, 
located 0.2° away from the phase center) not subtracted in the above step were imaged 
simultaneously from every time sample in the data to confirm the amplitude stability of 
non-variable sources. The rms noise in the maps are ~34 mJy.beam -1 in each polarization. 
The equivalent Stokes I noise ~25 mJy.beam -1 is approximately a factor of 2 lower than in 
Hyman et al. 2006. We attribute the significant reduction in noise due to better removal of 
extended sources before imaging GCRT coupled with more aggressive inner uv cutoffs (3.3 
k-Lambda vs. 2 k-Lambda employed by Hyman et al. 2006). 
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3. Results 

Light curves of the unresolved source, the GCRT, during its active state on September 
28, 2003 along with the unrelated compact field source G358. 917+0. 073 using their peak 
flux densities are shown in Fig. 2 (top panel) for both the polarizations. 

600 
^ 500 



I 400 



E 300 

to 

a 200 

55 100 

80 

60 
40 
20 


-20 



> 

CD 



O 

co -40 
-60 
-80 
-100 



CD 
O 
CO 



0.5 



-0.5 



G358.91 7+0.91 7 (RR/2.5) 
G358.91 7+0.91 7 (LL/2.5) ---a-- 
GCRT RR i— i— r 
I GCRT LL 



50 100 150 200 

Time (sec); start 28 Sept 2003 1 1 :45:00 (I AT) 



Fig. 2. — Light curves of the GCRT and an unrelated compact field source G358. 917+0. 073 
scaled down by 2.5 for both RR and LL (top). The arrows are 3a upper limits for non- 
detections. Other 2 panels show Stokes V at the same timestamps as above (middle) and 
the ratio of Stokes V to I (bottom). 



Due to a correction in data timestamps produced by upgraded observatory software, 
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the time used in this paper lags by 7.4 sec relative to iHyman et al.l ( 120061 ). The peak flux 
density of both the sources are smaller for the first data point, which is erroneously reduced 
by the data acquisition software at the start of data acquisition scans. The plot shows that 
while the flux density for both the polarizations of the unrelated field source are equal and 
steady within the error-bars, the flux density of the GCRT varies systematically with time. 
At the start of the scan, the RR flux density is higher than the LL, but it drops almost 
exponentially to zero after about 200 s econds and after wh ich the GCRT could not be 



identified. Due to lower rms noise than 



Hyman et al 



(120061 ) . we have been able to detect 



the source for ~1.5 times longer duration than reported earlier. 



Hyman et al. 



( 120061 ) searched for circular polarization from this source in a map 
made using data averaged over the time scale of emission and derived an upper limit of 
15%. However, because the difference between the two polarizations changes sign within 
the averaging time, their ability to constrain stokes V was smeared out. Adopting their 
approach during our re-analysis yields a stokes V of 12%, consistent wi th their reported 



20091 ) yields an upper 



limit. More recent reanalysis of its 2002 outburst (jSpreeuw et al. 
limit of 8% on the fraction of circular polarization averaged over the emission timescale. To 
study the variation of the circular polarization, we have also made maps of the GCRT in 
Stokes V, and a light curve made from these maps is shown in Fig. 2 (middle panel). In 
Stokes V, the rms noise in the map is ~17 mJy.beam -1 , quite close to the expected thermal 
noise. 

The figure shows that the sign of Stokes V reverses after about 100 sec from the start of the 
scan and within the data averaging time of ~17 sec. If the fastest variation in the source 
properties approach the velocity of light, the emission region would be smaller than 8 R . 
Fig. 2 (bottom panel) shows the ratio of Stokes V to I obtained by dividing the Stokes V 
maps by the Stokes I maps of the same timestamps using COMB in Aips. It shows the 
fractional circular polarization is initially low, but increases with time. After 100 seconds 
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from the start of the scan it reverses sign. The fraction of circular polarization approaches 
~100% about 150 seconds from the start of the scan. Fig. 3 shows a gray scale map of the 
difference of Stokes V images made from data integrated between 35 to 103 and 103 to 155 
seconds from the start of scan, respectively, during which the sign of stoves V reversed. This 
result (>6 a) represents our highest significance detection of Stokes-V from the GCRT. 



GREY: G358.7+0 VPOL 333.125 MHZ GCRT03V.DIFF.VDIFF.1 

-50 50 100 




RIGHT ASCENSION (J2000) 
Grey scale flux ranges -74.6 121.5 MilliJY/BEAM 

Fig. 3. — Stokes V difference image made 
from data integrated between 35 to 103 and 
103 to 155 seconds respectively from the 
start of scan and bridging the polarization 
reversal. Rms noise 19.2 mJy.Beam -1 . The 
location of the GCRT is at the center. 



Due to large closure errors on the calibrators with lower side band data from the 
correlator, the bandwidth of the present analysis was limited to 16 MHz. Th e measured 



Stokes I spectral index (S oc v a ) is —4 ± 3, which is better constrained than in 



Hyman et al. 



( 120061 ) (-4 ±5). We have also reanalyzed the September 2002 and March 2004 GCRT 
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detections. However, we could not detect any significant circul ar pol arization from those 2 
epochs, and our upper limits remain as in 
respectively. 



Hvman et al. 


(2005 


) and 


Hvman et al. 


(2007) 



4. Discussions 



4.1. Emission mechanism 



Fr om the measu r ed flu x density and the decay time scale (~ 2 min) of the 2002 



bursts, 



Hvman et al 



14 pc ( jSpreeuw et al. 



i 2005) placed an upper limit of 70 pc, which has been revised to 
20091 ) for incoherent emission. Adopting the ~ 17 sec time scale 
for polarization reversal derived from our new analysis of the detected 2003 burst as the 
fastest variation yields an upper limit of 8 R Q on its size. Assuming the above size of 
emission region in 2002 emission lowers the distance limit to ~ 10 pc. Given a very low 
probability of encountering an object within 1 pc and taking this as a distance lower limit 
still implies a brightness temperature >10 10 K. Such a high brightness temperature and 
circular polarization fraction reaching ~100% rules out ther mal emission and incoh erent 



19821), and 



non-thermal emission like synchrotron and gyro-synchrotron (IMelrose fc Dulkl 
indicates a coherent emission mechanism. Coherent radiation is known to arise from 3 
different processes, (i) pulsar emission, (ii) electron cyclotron maser emission and (iii) 
plasma emission, which are discussed below. 

Pulsar emission could produce highly circularly polarized emission within a single 
pulse. However, the integration time (16.8 sec) used in the data is longer than the pulse 
period of any known pulsar, and circular polarization fraction of pulsars goes down when 
averaged over a full pulse and is typically found to be significantly less than 100% near an 
observing frequency of 325 MHz (D. Mitra, private communication). This argues against 
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any pulsar-based model as an explanation for the observed properties of the GCRT. 
Electron cyclotron maser emission occurs at the electron gyr o-frequency and its harmonics 



19851). If the GCRT 



and could also produce highly circularly polarized emission ( iDulkl 
emission at 330 MHz was due to cyclotron emission at its fundamental frequency, the 
required magnetic field is ~120 Gauss and ~2 lower at the second harmonic. These 
are typical values of magnetic field in a stellar corona. The relative bandwidth of 



emission is Avjjs re Vn/c 2 



(Ivan den Oord &: de Bruyn 



or e 



ectrons with speed v , and is typically ~ 0.01 — 0.1 
19941 ). We note, however, that the individual spikes of emission 
could have durations much smaller than our integration time and could have occured from 
different physical locations. A significant variation of magnetic fields within these spikes of 
emission would cause a variation of emitting frequencies causing a br oadening of observed 



bandwidth. This could account for >30 MHz bandwidt 



i observed in 



The very steep spectral index reported in iHyman et al 



Hyman et al. 



(I2007f ). 



(120071 ) could also be explained if 



the peak of emission at that time occured at <310 MHz. 



Plasma emission occurs at the plasma frequency or one of its lower order harmonics. 
It involves plasma instabilities and subsequent conversion of a fraction of its energy to 
electro magnetic radiation. Different me chanisms of generating plasma instabilities are 



known (Ivan den Oord fc de Bruyn 



19941 ). Common characteristics of these are narrow 



band emission with a high level (approaching 100%) of circular polarization and high 
brightness temperature at meter wavelengths. The typical bandwidth of plasma emission is 
Av/v ~ 0.01. However, as described above, there could be several spots of emission slightly 
differing in fr e quenc y, which would also explain a bandwidth of >32 MHz as reported in 



Hyman et al 



(120071 ). We note that plasma frequency of 330 MHz corresponds to an ion 
density of ~10 9 ions. cm -3 , which is less likely but cannot be ruled out in the corona of a 



brown dwarf. However, untill the magnetic field strengths or ion density could be measured 
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near the object, any electron cyclotron or plasma emission mechanism attributable to the 
object cannot as yet be firmly established. 



4.2. Distance to the GCRT and its classification 



If the upper limit to the brightness temperature (Tb) and the linear size (L) of the 
emission region is known, an upper limit to its distance (D) can be calculated. For electron 
cyclotron maser emission, it is found that Tb could reach 10 20 K as in auroral kilometric 



(Me 



1999 



Nik 



radiation (lErgun et al. 



2000). This upper limit has also been estimated t heoretically 



19941 ) . For plasma emission, T B is expected to be less than 10 16 K (IStepanov et al. 



However, 


'or a 


(Kuijpers 


1989) 



or a certain type of plasma emission from double layers, Tb could 



19891 ) . Therefore, we could easily consider the upper limit on its Tb 
to be 10 20 K. With L <8 R (Sect. 4), the upper limit on its distance is ~ 100 kpc, which 
indicates that the emission does not originate from another galaxy. However, this limit can 
be much improved by considering known sources of circular polarisation in the Galaxy, 
which are not ruled out already. 

Most of the known low frequency emitting objects are either synchrotron or thermal 
sources, from which the circular polarization fraction is less than a few percent. However, 
flare stars in the Galaxy could emit strong high circularly polarized emission. These objects 
are typically dwarfs stars of class G to M (mostly M). Som e of them c o uld b e of size ~0.15 



R (e.g., UV Cet with spectral type M6). As described in iGershbergl ( 120051 ). more than 
50 observations have been carried out in radio bands towards them from the early days of 
radio astronomy in 1958 to 2001 and cover a frequency range from 20 MHz to several GHz. 
The flux densities recorded for these flares at cm wavelengths are typically in millijanskys, 
and the highest recorded was only a fraction of a Jy. However, at lower frequencies, 
particularly at decameter wavelengths, flux densities of some of the flares recorded are 
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100 Jy (lAbdul-Aziz et al.lll995h to ~ 10 4 Jy fjSleel Il963h . The highest lu minosity of any 



flare s tar(s) in the metre wave band is ~230 Jy at 136 MHz as recorded by 



Slee fc Higgins 



( 119691 ) towards the Orion nebula located ~400 pc away If the luminosity of the GCRT in 



the 2002 outburst is the same, it would be located at a distance of ~4 kpc. 

Earlier observations of flare stars showed wide variation in the circular polarization 
fraction with some of the bursts being unpolarized while others h aving circular polarization 



fraction reaching up to 100% (see e.g., 



Nelson et al. 



stars. 



Abada-Simon et al. 



Lang et al. 



(119941 ) ) . At meter wavelengths, 



( 1979) obs erved 40 to 60% circularly polarized emission in 8 cases from flare 



( 119831 ) reported observations of AD Leo with circular polarization fraction 
of about 15%, and the emission properties were explained as electron cyclotron emission. 
Reversal in the sens e of circular pola rization was also reported from the above flare star 



( jJackson et al. 



19891 ). As discussed in 



DulM (jl985h and 



Melrose! ( 1l980l ). mode coupling could 



prevent attainment of a high polarization, which could operate for the GCRT when the 
circular polarization fraction is small and also explain our observed polarization reversal in 
the data described above. 

Measured Rotation period for stars show a decrease in r otation period to wards lower 



mass (0.1 M ) objects reaching ~0.1 day (e.g. in Pleaide s) (llrwin et al. 



th ese ultracool dwarfs p roduce flares at radio frequencies ( IBurgasser &: Putman 



2008). A few of 



in 



Hallinan et al. 



20051). As 



( 120071 ) . within one rotation of the object, there could be two pulses of 
emission. If this is the case during the GCRT 2002 outbursts, its rotation period is ~0.1 
day, consist ent with the above finding. Hence, an ultracool dwarf could be progenitor of 



the GCRT. 



Kaplan et al. 



( 120081 ) have considered the object 'C as a possible progenitor of 
the GCRT. However, emission from an ultracool dwarf (L4.5V) do not produce a good fit 
to the observed multi-band photometry of the object C. We find the goodness of the fit do 



not change significantly (rms error ~0.3 magnitude) from 



Kaplan et al. 



(120081 ) results when 
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the suggested K7V star is replaced by a ~0.1 M star ~0.1 Gyr old (iBaraffe et all 119981 ). 
This could be a young mid to late M type of star with Ay ~3.5, indicating it to be at a 
distance of ~4 kpc. Such a distance matches the 10 20 K cyclotron brightness temperature 
limit (radius ~0 .15 R^). However, considering the uncertainty (~0.1 magnitude) in the 



IR magnitudes ( iKaplan et al 



20081 ) and in the theoretical model, the quality of the fit is 



marginal and object 'C is considered only as a candidate progenitor of the GCRT. We, 



howe ver, note that a hithertho undetected ultracool dwarf located >200 pc ( IKaplan et al. 



20081 ) could also be the progenitor of the GCRT. 



5. Conclusions 

We have detected time-varying, highly circularly polarised emission from the GCRT 
J1745-3009 based on a re-analysis of its outburst on 28th September 2003. The percentage of 
circular polarization approaching ~100% commensurate with high brightness temperature 
(>10 10 K for a source >1 pc away) significantly strengthens the original suggestion of a 



coherent emission scenario ( Hyman et al. 



20051 ). but its properties are inconsistent with 



a pulsar origin. The time scale of change in sign of the circular polarization allows us to 
estimate an upper limit of 8 R on its size. Based on its measured emission characteristics, 
either electron cyclotron maser or plasma emission processes could be the most plausible 
origin scenarios. 

Ascribing the 77 minute periodicity observed in its 2002 outbursts to half the rotation 
period of its progenitor, observed radio emission of the GCRT could be explained as 
outbursts from a highly subsolar flare star. Its distance is estimated to be <4 kpc. 
Independent verification, would however, be required to confirm its progenitor. 
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